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ABSTRACT 



In our recent paper (lYadav et al.ll2007l ) we described a fast cubic (bispectrum) 



estimator of the amplitude of primordial non-Gaussianity of local type, /nl, from 
a combined analysis of the Cosmic Microwave Background (CMB) temperature 
and E-polarization observations. In this paper we generalize the estimator to 
deal with a partial sky coverage as well as inhomogeneous noise. Our generalized 
estimator is still computationally efficient, scaling as 0{N^l'^) compared to the 
0{Npl^) scaling of the brute force bispectrum calculation for sky maps with 
Npix pixels. Upcoming CMB experiments are expected to yield high-sensitivity 
temperature and E-polarization data. Our generalized estimator will allow us to 
optimally utilize the combined CMB temperature and E-polarization information 
from these realistic experiments, and to constrain primordial non-Gaussianity. 

Subject headings: cosmic microwave background, early universe, inflation, non- 
Gaussianity. 
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1. Introduction 



Non-Gaussianity fro m the simplest inflation models, that are based on a slowly rolling 



scalar fleld, is very small flSalopek fc Bondl 11991 Il99ll : iFalk et al.lll993l : iGangui et all 11994 
Acquaviva et al. 2003 : Maldacenal 2003 ): however, a very large class of more general models, 



e.g., models with multiple scalar flelds, features in inflation potential, non-adiabatic fluctua- 
tions, non-canonical kinetic terms, deviations from the Bunch-Davi es vacuum, among o thers, 
predict substantially higher level of primordial non-Gaussianity (IBartolo et al.l l2004j . for a 
review and references therein). 

Primordial non-Gaussianity can be described in terms of the 3-point correlation function 
of Bardeen's curvature perturbations, $(fc), in Fourier space: 

($(ki)(k2)(k3)) = i2nf5%k^ + k2 + ks)Fih, h, h). (1) 

Depending on the shape of the 3-p oint function, i.e., F {k^, ^2, k^)-, non-Gaussianity can be 
broadly classifled into two classes (IBabich et al.ll2004l ). First, the local, "squeezed," non- 
Gaussianity where F{ki, k2, k^) is large for the conflgurations in which ki <C k2, k^. Second, 
the non-local, "equilateral," non-Gaussianity where F{ki, k2, k^) is large for the conflguration 
when A;i ~ A;2 ~ ^3- 

The lo cal form arises from a non-l i near relation betwe en inflaton and cur vature per 



turbations (ISalopek fc Bondl Il990l . Il99ll: iGangui et a. 



2003), or the New Ekpyrotic models (jKoyama et al. 200?! : Buchbinder et al. 20071 ). The 



1994 ) , curvaton models (ILvth et al. 



equilat eral form arises from no n-canonical kinetic terms s uch as the Dirac-Born-Infe ld (DBl) 
action (lAlishahiha et al.ll2004j ). the ghost condensation (lArkani-Hamed et al.l 120041). or any 



other single-field models i n which the scalar field acquires a low speed of sound (jChen et al. 



20071 : ICheung et al.l 120071 ). While we focus on the local form in this paper, it is straightfor- 
ward to repeat our analysis for the equilateral form. 



The local form of non-Gaussianity may be pa rametrized in real space as (IGangui et al. 



1994J : I Verde et al.ll2000l : iKomatsu fc Spergelll200lh : 

$(r) = $L(r) + /^z.(«fi(r)-(«fi(r))) (2) 

where /atl characterizes the amplitude of primordial non-Gaussianity. Different inflationary 
models predict different amounts of /atl, starting from 0(1) to f'NL ~ 100, beyond which val- 
ues have been excluded by the Cosmic Microwave Bac kground (CMB) bispectrum of WMAP 
tempe r ature data, —36 < J nt, < 100, at the 2a level ( iKomatsu et al.ll2003l : ICreminelli et al. 
2006bl : ISpergel et al.l 120061 ). 



So far all the constraints on primordial non-Gaussianity use only temperature infor- 
mation of the CMB. By also having the E-polarization information together with CMB 
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temperature information, one can improve the sensitivity to the prim ordial fluctuations 
teabich fc Zaldarriagalbooi : IVadav fc Wandeltl boosl : IVadav et al.lboOTi ). Although the ex- 



periments have already started characterizing E-polariza tion anisotropies (iKovac et al.ll2002 



Kogut et al.ll2003l : iPage et al.l 120071 : iMontroy et al.ll2006l ). the errors are large in comparison 
to temperature anisotropy. The upcoming experiments such as Planck satellite will charac- 
terize E-polarization anisotropy to a higher accuracy. It is very timely to develop the tools 
which can optimally utilize the combined CMB temperature and E-polarization information 
to constrain models of the early universe. 

Throughout this paper we use the standard Lambda CDM cosmology with the following 
cosmological parameters: fib = 0.042, Qcdm = 0.239, Ql = 0.719, h = 0.73, n = 1, and 
r = 0.09. For all of our simulations we used HEALPix maps with Npi^ !^ 3 x 10® pixels. 



1.1. Generalized Bispectrum Estimator of Primordial Non-Gaussianity 



In our recent paper (lYadav et al.ll2007l ) we described a fast cubic (bispectrum) estimator 
of /nl, using a combined analysis of the temperature and E-polarization observations. The 
estimator was optimal for homogeneous noise, where optimality was defined by saturation 
of the Fisher matrix bound. 

In this paper we generalize our previous estimator of /nl to deal more optimally with a 
partial sky cover age and the inho r nogene ous noise. The generalization is done in an analo- 
gous wa y to how ICreminelli et al.l (j2006al ) generalized the te mperature only estimator devel- 
oped by iKomatsu et al.l (120051 ): however, the final result of ICreminelli et al.l (l2006al ) (their 
Eq. (30)) is off by a factor of two, which results in the error in /tvl that is much larger than 
the Fisher matrix prediction, as we shall show below. 

The fast bispectrum estimator of Jnl from the c ombined CMB tem perature and E- 
polarization data can be written as /nl 

1 



prim 



fsky 



r dr 



where ( lYadav et al.ll2007h 

(Pn B'^{h,r)A(h,r), 



(3) 



E E 

ijkpqr 2<£i<£2<^3 



A, 



Tjpqr,prim / ^-l\ip / ^-l\jq / ^-l\kr n 

'^iii2i3 )^2y^ hi^. 



ijk,prim 



(4) 
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(5) 
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A{n,r) = J2J2iC~'^'''"' 

ip Im 



(6) 



TT 



edkP^{k)g\{k)ji{kr) 



(7) 



-i j kHk gmUkr) 



and fsky is a fraction of the sky observed. Indices i,j,k,p,q and r can either be T or i?. 
Here, Ag-^g^e^ is 1 when ii £2 7^ £3 , 6 when ^i = £2 = ^3, and 2 otherwise, B^^^l'^^^ is the 
theoretical bispectrum for f^i = 1 (jYadav et ahl 120071 ). P<s>{k) is the power spectrum of the 
primordial curvature perturbations, and g}{r) is the radiation transfer function of adiabatic 
perturbations. 

It has been shown th at the above menti oned estimator is optimal for the full sky coverage 
and homogeneous noise (jYadav et al.l 120071 ). To be able to deal with the realistic data, the 
estimator has to be able to deal with the inhomogeneous noise and foreground masks. 

The estimator can be generalized to deal with a partial sky coverage and the inhomoge- 



neous noise by adding a linear t erm to .S*; 



prim. '■ Snrim. 



linear 
prim ■ 



For the temperature 

only case, this has been done in ICreminelli et al.l (j2006al ). Following the same argument, we 
find that the linear term for the combined analysis of CMB temperature and polarization 
data is given by 



5': 



linear 
prim 



fsky 



r'^dr / d^n {2B{n, r) {Asim{n, r)Bsim{n, r))MC + -4(n, r) (5^j„(n, r))Mc} , (9) 



where Asim{n,r) and Bsim{'>T',r) are the A and B maps generated from Monte Carlo sim- 
ulations that contain signal and noise, and (..) denotes the average over the Monte Carlo 
simulations. 

The generalized estimator is given by 



NL 



q I qlinear 

^prim I '-'prim 

N ' 



(10) 



which is the main result of this paper. Note that {Sp^^"^) mc = —{Sprim) mc, and this relation 
also holds for the equilateral shape. Therefore, it is straightforward to find the generalized 
estimator for the equilateral shape: first, find the cubic estimator of the equilateral shape. 
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Sequilaterah and take the Monte Carlo average, {S equilateral) MC- Let us suppose that SequUateral 

contains terms in the form of ABC, where A, B, and C are some filtered maps. Use the 
Wick's theorem to re- write the average of a cubic product as {ABC)mc = {A)mc{BC)mc + 
{B)mc{AC)mc + {C) Mc{-^B) MC- Finally, remove the MC average from single maps, and 
replace maps in the product with the simulated maps: {A)mc{BC)mc + {B)mc{AC)mc + 
{C)mc{AB)mc A{BsimCsim) MC + B {AsimC sim) MC + C {AsimB sim) MC ■ This Operation 
gives the correct expression for the linear term, both for the local form and the equilateral 
form. 

One can find the estimator of f^L from the temperature data only by setting i = 
j = k = p = q = r = T. We hav e compared our f o rmula in the temperature-only limit 
with the original formula derived by ICreminelli et al.l (l2006a( ) (their Eq. (30)), and found a 

discrepancy. To see the discrepancy, let u s re-write the est i mator as: f^i 
Our formula gives x = 1, while Eq. (30) of ICreminelli et al.l (j2006al ) gives x = 0.5. 



-iN 



To make sure that our normalization gives the minimum variance estimator, we have 
done Monte Carlo simulations with varying x. We find tha t x = 1 minimiz es the variance (as 



shown in Fig. [T]). We conclude that the analysis given in ICreminelli et al.l (l2006al ) resulted 
in the larger-than-expected uncertainty in /^vl because of this error in their normalization 
of the linear term. 

The main contribution to the linear term comes from the inhomogeneous noise and 
sky cut. For the temperature only case, most of the contribution to the linear term comes 
from the inhomogeneous noise, and the partial sky coverage does not contribute much to 
the linear term. This is because the sky-cut induces a monopole contribution outside the 
mask. In the analysis one subtracts the monopole from outside th e mask before measurin g 



Sprimi which makes the linear contribution from the mask small (ICreminelli et al.ll2006al ). 



For a combined analysis of the temperature and polarization maps, however, the linear 
term does get a significant contribution from a partial sky coverage (see the right panel of 
Fig. [2]). Subtraction of the monopole outside of the mask is of no help for polarization, as 
the monopole does not exist in the polarization maps by definition. (The lowest relevant 
multipole for polarization is / = 2.) 

The estimator is still computationally efficient, taking only N'^i'^ (times the r sampling, 
which is of order 100) operations in comparison to the full bispectrum calculation which takes 
^pix operations. Here Npi^ refers to the total number of pixels. For Planck, Npi^ ~ 5 x 10^, 



^Equation (30) in ICreminelli et al.1 (|2006al) is off by a factor of 2. Since w e used X^^i^j^a ~ 
6 X]fi<^2<^3 a7T1~' compare our normalization factor x with ICreminelli et al.l (|2006a| ) one needs to divide 
the normalization in ICreminelli et al.l (|2006al) by a factor of 6. 
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100 





Normalization factor (x) 



Fig. 1. — Testing normalization of the linear term in the estimator of Jml- The symbols show 
the standard deviation of Jni derived from the the Monte Carlo simulations using the esti- 
mator for a given normalization , x. The horizontal line s hows the Fisher matrix prediction. 
Our formula gives x = 1, while ICreminelli et al.l (j2006al ) give x = 0.5 (their Eq. (30)). We 
have used simulated polarized Gaussian CMB maps with the Planck inhomogeneous noise 
as well as the WMAP KpO and P06 mask for temperature and polarization, respectively. 



and so the full bispectrum analysis is not feasible while our analysis is. 



2. Results 

In the left panel of figure [21 we show the variance of fjyi using the estimator (with and 
without the linear term) for the Gaussian CMB simulations in the presence of inhomogeneous 
noise and partial sky coverage. For this analysis we use the noise pro perties that are expecte d 
for the Planck satellite, assuming the cycloidal scanning strategy (IDupac fc Tauberl l2005l ) . 
Inhomogeneous nature of the noise is depicted in the lower map of figure H] where we show 
the number of observations (Nobs) for the different pixels in the sky. As for the foreground 
masks, we use WMAP KpO intensity mask and P06 polarization mask. We find that, with 
the inclusion of the linear term, the variance reduces by more than a factor of 5. The 
linear term greatly reduces the variance approaching the Fisher matrix bound. However, the 
estimator is close to, but not exactly the same as the Fisher variance prediction in the noise 
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Fig. 2. — Optimality of the generalized estimator. The sohd hnes show the Fisher matrix 
prediction for the standard deviation of f^L, the triangles show the standard deviation 
derived from the Monte Carlo simulations using the estimator without the linear term, and 
the stars show the standard deviation derived from the Monte Carlo simulations using the 
generalized estimator (i.e., with the linear term). Left panel: The uncertainty vs the 
maximum multipole that is used in the analysis, £max- The simulations contain the Gaussian 
CMB signal, inhomogeneous noise (which simulates the Planck satellite), WMAP KpO and 
P06 masks. Right panel: The uncertainty vs a fraction of the sky observed, fsky, for 
^max = 500. The simulations include the Gaussian CMB signal, and flat sky-cut (which is 
azimuthally symmetric in the Galactic coordinates), while they do not include instrumental 
noise. This figure therefore shows that the sky cut contributes significantly to the linear 
term of polarization. 



dominated regime. 



Nevertheless, we do not observe the increase of variance at higher Imax'- the variance 
be comes smaller as we inclu de m ore multipoles. This resu lt is in contradiction with the result 
of ICreminelli et al.l (l2006al ) and ICreminelli et al.l (l2006bl ) . We attribute this discrepancy to 
the error in the normalization of linear term in their formula. 

In the right panel of figure [2l we show the variance of /tvl again using Gaussian sim- 
ulations, but now in the presence of flat sky cut and in the absence of any noise. The 
purpose of the plot is to demonstrate (as pointed out in the previous section) that for the 
combined CMB temperature and polarization analysis, sky-cut does contribute significantly 
to the linear term. We find that the generalized estimator does a very good job in reducing 
the variance excess, and the simulated variance of /tvl does accurately saturate the Fisher 
matrix bound. 
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Noise 


Sky-cut 


(fNL) 


ninput 


^sim 


No 


flat cut, /sfcy = 0.8 


103.2 


100 


10.1 


Inhomogeneous 


WMAP KpO and P06 masks 


108.7 


100 


21.04 



Table 1: Unbiasedness of the generalized estimator. Non-Gaussian CMB maps with fjjj'jf"* = 
100 are used for i^ax = 500. The standard deviation of /atl, agim, was obtained using 
Gaussian simulations. 



<A(r)B(r)>„,- <B(r)B(r)>„ 




Fig. 3. — The {Asim{n,r)Bsira{n,r)) mc and {Bl^^{h^r)) mc maps in dimension-less units 
for a slice near the surface of last scattering. These maps are calculated from Monte Carlo 
simulations with the Gaussian signal, Planck inhomogeneous noise, and WMAP KpO and 
P06 masks. 



Can our estimator recover the correct f^L-, i-e., is our estimator unbiased? We have 
tested our estimator against simulated non-Gaussian CMB temperature and E-polarization 
maps. The non-Gaussia n CMB temperature and E-polarization maps were generated using 
the method described in iLiguori et al.l (120071 ). We find that our estimator is unbiased, i.e., 
we can recover the /tvl value which was used to generate the non-Gaussian CMB maps. The 
results for the unbiasedness of the estimator a re shown in Table [T] The analysis also shows 
the unbiasedness of the estimator described in lYadav et al.l (120071 ). 



Figures E] and H show the maps {Asim{n,r)Bsirn{n,r)) mc and {Bl-^{h,r)) mc, which 
appear in the linear term (Eq. [H]) of the estimator. These maps are calculated using 100 
Monte Carlo simulations of the data. Since the linear term contributes only in the presence 
of inhomogeneities, we also show these maps calculated with noise-only simulations (i.e. no 
signal). Notice how these maps correlate with the inhomogeneous noise (as shown in the 
lower map of figure Sj) . 
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<A(r)B(r)>„, 



<B{r)B(r)>„ 





i^^^^^^^^^^^^^^^^H 61374 N„^, 



Fig. 4 — The top panels show {Asim{n,r)Bsi,n{n,r)) mc and {Bl^^{h,r)) mc maps for the 
noise only analysis (i.e. no CMB signal or mask). The maps are in dimension-less units and 
are shown for a slice near the surface of last scattering. The Lower map shows the number 
of observations per pixel (Nobs) at the resolution of Npi^ = 12582912. 



Conclusion and discussion 



Upcoming CMB experiments will provide a wealth of information about the CMB po- 
larization anisotropies together with temperature anisotropies. The combined information 
from the CM B temperature and polarizat i on data improves the sen s itivity to primordi al non- 
Gaussianitv JBabich fc Zaldarriagal booi IVadav fc Wandeltl boosl : IVadav et al.l[2007[ ). The 
promise of learning about the early universe by constraining the amplitude of primordial 
non-Gaussianity is now well established. I n this paper we hav e generalized the bispectrum 
estimator of non-Gaussianity described in lYadav et al.l (120071 ). to deal with the inhomoge- 
neous nature of noise and incomplete sky coverage. 



The generalization from lYadav et al.l (120071 ) enables us to increase optimality of the 
estimator significantly, without compromising the computational efficiency of the estimator: 



the estimator is still computa tionally efficient, scaling as 0(A^p/5) compared to the 0{N^l 



r5/2x 



jif I II ' ' I ^ pix / \ jjti. 

Babich fc Zaldarriagall2004l ) calculation for sky maps with 



pix J 
'I pix 



pixels. For the Planck satellite, this translates into a speed-up by factors of millions, reducing 
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the required computing time from thousands of years to just hours and thus making /nl 
estimation feasible. The speed of our estimator allows us to study its statistical properties 
using Monte Carlo simulations. 

We have used Gaussian and non-Gaussian simulations to characterize the estimator. We 
have shown that the generalized fast estimator is able to deal with the partial sky coverage 
very well and in fact the variance of f^i saturates the Fisher matrix bound. In the presence 
of both the realistic noise and galactic ma sk, we find that the generalized estimator greatly 



reduces the variance in comparison to the lYadav et al.l (120071 ) estimator of non-Gaussianity 



using combined CMB temperature and polarization data. 

Since the estimator is able to deal with the partial sky coverage very effectively, the 
estimator can also be used to constrain primordial non-Gaussianity using the data from 
ground and balloon based CMB experi ments which observ e only a small fraction of the sky. 



The estimator also solves the problem (lYadav et al.l 120071 ) of non-trivial polarization mode 



coupling due to foreground masks. Earlier this issue was dealt with by removing the most 
contaminated i modes from the analysis (usually £ < 30). 

The naive approach of using galactic masks to deal with the polarization contamination 
is to be refined. Both temperature and polarization foregrounds are expected to produce 
non-Gaussian signals. Some sources of non-primordial non-Gaussianity are CMB lensing, 
point sources, and the Sunyaev Zel'dovich effect. Understanding the non-Gaussianity from 
the polarization foreground sources and refining the estimator to be able to deal with it will 
be the subject of our future work. 



Some of the results in this paper have been derived using t he CMBFAST package by 
Uros Seljak and Matias Zaldarriaga (jSeljak fc Zaldarriagal 119961 ) and the HEALPix pack- 
age (iGorski et al.l 120051 ). This work was partially supported by the National Center for 
Supercomputing Applications under TG-MCA04T015 and by University of Illinois. We 
also utilized the Teragrid Cluster (www.teragrid.org) at NCSA. BDW acknowledges the 
Friedrich Wilhelm Bessel research award by the Alexander von Humboldt foundation. BDW 
and APSY also thank the Max Planck Institute for Astrophysics for hospitality. BDW and 
APSY are supported in part by NSF grant numbers AST 0507676 and 0708849, NASA/JPL 
subcontract no. 1236748. EK acknowledges support from the Alfred P. Sloan Foundation. 
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